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ABSTRACT 



Context. Weak gravitational lensing is a powerful probe of large-scale structure and cosmology. Most commonly, second-order corre- 
lations of observed galaxy ellipticities are expressed as a projection of the matter power spectrum, corresponding to the lowest-order 
approximation between the projected and the three-dimensional power spectrum. 

Aims. The dominant lensing-only contribution beyond the zero-order approximation is the reduced shear, which takes into account 
not only lensing-induced distortions but also the isotropic magnification of galaxy images. This involves an integral over the matter 
bispectrum. We provide a fast and general way to calculate this correction term. 

Methods. Using a model for the matter bispectrum, we fit elementary functions to the reduced-shear contribution and its derivatives 
with respect to cosmological parameters. The dependence on cosmology is encompassed in a Taylor-expansion around a fiducial 
model. 

Results. Within a region in parameter space comprising the WMAP7 68% error ellipsoid, the total reduced-shear power spectrum 
(shear plus fitted reduced-shear correction) is accurate to 1% (2%) for t < 10 4 (€ < 2 X 10 5 ). This corresponds to a factor of four 
reduction in the bias compared to the case where no correction is used. This precision is necessary to match the accuracy of current 
non-linear power spectrum predictions from numerical simulations. 

Key words, cosmology - gravitational lensing - large-scale structure of the Universe 
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1. Introduction 

Cosmic shear surveys exploit the distortions of distant galaxy 
images induced by the matter structures in the Universe on large 
scales. Weak cosmological lensing has become an important tool 
to measure cosmological parameters. The current precision in 
the measured shapes of high-redshift galaxies is at the few per- 
cent level providing important constraints on the matter density 
Q m and the power-spectrum normalisation cr 8 . Future surveys 
aim to achieve sub-percent level accuracy in measuring dark- 
energy properties or deviations from general relativity. 

To reach this precision for cosmological parameters, efforts 
not only on the observational but also on the theoretical side have 
to be made. To predict the shear power spectrum from theoret- 
ical models, several approximations are usually made to facili- 
tate calculations. When comparing and fitting models to obser- 
vations, these simplifications can significantly bias the resulting 
cosmological parameters. One of these approximations is to ne- 
glect the correction due to the so-called reduced shear. 

The observed ellipticity of galaxies is an estimator of the 
reduced shear 



r 



l-K 



(1) 



where y is the shear and k is the convergence. Here and in the fol- 
lowing, both shear and reduced shear are represented as complex 



quantities, y - j\ + iyi, g — gi + igi- Since in weak lensing the 
convergence is small, \k\ <s: 1, one usually simplifies the above 
relation to linear order, g ~ y. 

The effect o f reduc e d she ar was examined using /V-body 
simulations by IWhitd (20051) an d com pared with theoreti- 
cal predictions in Dodelson et al.l (l2006|). Corrections to the 
lensing skewn ess were disc ussed in Schneid er et a D (119981) 
and iDodelson & Zhang d2005l) . An extensive study of higher- 
order corrections to the weak-lensing power spectrum, includ- 
ing reduced shear, was presented by iRrause & Hiratal ([2009). 
Compared with corrections due to light-path integration (relax- 
ation of the Born approximation) a nd lens-lens coupling (e.g., 
ISeitzl [19931 iBernardeau et aD 1 19971) . reduced-shear corrections 
represent the dominant contribution to the weak-lensing power 
spectrum. The lowest-order correction term, which is of order 
(9(0 3 ) in the gravitational potential <p, reaches 10% of the to- 
tal lensing spectrum. This is the term that we consider in this 
work. The next-orde r correction at Q(<ft 4 ) is one magnitude be- 
low the 0(0 3 )-term dKrause & Hiratal l2009h . Ignoring reduced 
shear will cause a bias in cosmological parameters, which is 
comparable to the statist ical error f rom future surveys likeDE^l 
LSSlQ, or EuclicQ ( Shapiro 2009). For the weak-lensing results 
of the COSMOS survejQ an estimate from numerical simu- 
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lation s yielded a 1% underestimation of a% dSchrabbac k et al.l 
I2010h . 

Another correction to weak lensing from a magnitude- 
limited galaxy sample is the so-called magnification bias. 
iHamanal (12001) found its contribution to the lensing power spec- 
trum to be small. The corresponding third-order correction terms 
we re also calcul a ted an d compared to numerical simulations 
by iMenard et al.l (120031) . However, if not corrected for, future 
surveys wil l yield significantly b iased cosmological parameter 
constraints (Schmidt etal. 2009a). The lowest-order correction 
term has the same functional form as the reduced-shear one, and 
therefore both can be dealt with simultaneously. Higher-order 
terms cont ribute by smaller than o ne percent to the lensing power 
spectrum dKrause & Hiratall2009l) . 

There are additional important sources of contamination to 
the weak lensing observables, some of which may be larger than 
the reduced-shear corrections. These involve galaxy clustering, 
the intrinsic alignment of galaxies with each other or with the 
surrounding dark matter haloes, and the modelling of the non- 
linear dark-matter power spectrum and baryonic processes. The 
modelling of these effects requires detailed understanding of 
baryonic physics, galaxy formation, the interplay between galax- 
ies and dark matter, and massive numerical simulations in the 
case of the non-linear power spectrum. In contrast, taking into 
account the reduced shear is straightforward and can be done 
with sufficient accuracy to render it a sub-dominant contamina- 
tion. 

In the next section, we review the expressions for the lowest- 
order reduced-shear term, and introduce a fitting scheme to ef- 
ficiently model this term. In Sect. [3] we present our results, 
demonstrating the accuracy of the fits. We conclude this paper 
in Sect.0 



2. Reduced shear 

2.1. Lowest-order expansion 

For small values of the convergence, |a:| < 1, the reduced shear 
(Q} is expanded as 



1 -k 



= y[l + k + /c + . . .]. 



(2) 



To calculate the power spectrum of the reduced shear g, we apply 
a Fourier transform to Eq. (f2| to second order 



In our notation, the convolution of y and k is 



(9 



r d 2 r 
J W 



y(t')Ul-t'). 



(4) 



We assume that the reduced shear is measured in redshift bins. 
With P™, we denote the reduced-shear (cross-)power spectrum 
for two redshift bins with indices m and n. In general, the shear 
field can be decomposed into an E-mode (gradient) and a B- 
mode (curl). However, as is the case for cosmic shear at low- 
est order, the reduced-shear correction of order 0( d> 3 ) does not 
have a B-mode component dSchneider et alj|2002l) . and in the 
following we focus only on the E-mode. In the absence of a 
B-mode, the reduced -shear power spectrum is defined to be 
dSchneider et alj|2002l) 



<rW*(0> = (2n) 2 6 D (e - t')P'""({), 



(5) 



where g" 1 denotes the reduced shear measured from galaxies in 
the m th redshift bin, the asterisk (*) denotes complex conjuga- 
tion, and 5d is the Dirac delta-function. 

The reduced-shear power spectrum is the sum of the conver- 
gence power spectrum (which equals the shear power spectrum) 
and a correction given by 



(6) 

(2tt) 2 5 d {{ - t')P m ":\() = (T(t){y n * *T(0> + (m ~ n). (7) 



Using the relation between shear and convergence, y(() = 
exp(2i/?)£(f ), where B is the polar angle of the wave vector ( , and 
that the convergence is a real field, k*(() = k(-(), the third-order 
correlator in Eq. ^} can be written in terms of the convergence 
bispectrum B K , which is defined by the equation 

(FXA W 2 )F(/3)> = (2n) 2 6 D (e i +e 2 + h)BT'\txJi). (8) 
With this, we derive the correction term 



f TTTi^' i B T"(fo, -f) + BTih, -(')} . (9) 
J \2ny 



Without loss of generality, we have set B — and introduced 
(q = ({, 0). In the integrand, B' is the polar angle of t' . 

The convergence bispectrum can be expressed in terms of the 
three -dimensional matter bispec trum Bg using Limber's equa- 
tion dLimbetlll953t lKaiserHl992h providing 



bT({ 



,h) = f 

•^lir 



da F nmp (a) B s 



f K [x(a)Y f K lx(a)] 



;*(«); 



F mnp (a) = 



d X G m (a)G n (a)G p (a) 
da f K [x(a)] 



(10) 



The lower limit of the integral over the scale factor a corresponds 
to the limiting redshift, an m = 1/(1 + Zi; m ), x i s me comoving 
distance, and fx is the comoving angular diameter distance. The 
lensing efficiency G m for the m redshift bin is 

G (a) = x — dap m (a) 777737; . (H) 

2 \ c I a J fl| . m fKlx(a')] 

which includes an integral over the probability distribution p m 
of the number density of source galaxies in the m th bin. 



(3) 2.2. Expansion around a fiducial model 



In the expressions needed to calculate the first-order reduced- 
shear power spectrum in Eqs. d9l[T0b. we can interchange the l- 
and the a-integration, and replace the latter by a finite sum. Then 



F" m (a) = i [F mm (a k ) + F""""(a k )] ; 

Q((,a) = 2 [ £L COS (2B')bJ /o 
J (2n) 2 



(12) 



Evaluating the ^'-integral is not a difficult but very time- 
consuming numerical problem. It is significantly slower than the 
calculation of the matter power spectrum necessary for the con- 
vergence power spectrum. This makes reduced-shear corrections 
the limiting factor in the parameter analysis of a weak-lensing 
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survey, for example using Monte Carlo sampling. In the follow- 
ing, we present fitting formulae that provide an accurate approx- 
imation for this term within a region in parameter space around 
a fiducial cosmological model. The ^-independent term F mnn (d) 
is easily computed and can be tabulated. Since it depends on the 
source redshift distribution, it has to be determined for each sur- 
vey. 

We perform a Taylor-expansion to first order P g as a func- 
tion of a cosmological parameter vector p = (pi, pi, ■ ■ ■ Pm) 
around a fiducial cosmological model with parameter po 

P^T^P) - P m T«,Po) + (P - Po) V P P m " g m (Ap)\ p __ po . (13) 

Inserting Eq. (fT2l) . the derivatives with respect to cosmological 
parameters are given by 



dP^" g "\(, P ) 

dp a 



k=\ 



= AaY J \F" m (.a k ,p~) 



dQ({,a k ,p) 

dp a 



dF"'"(a k ,p) 1 
+ — A±Pl Q{(ak A (14) 

Op a ) 

In Sect. [3] we consider the cosmological parameters p = 
(Q m , Q(j e , w, Q.b, ft, erg, n s ) (see also TableQ]). 



2.3. Magnification and size bias 

A typical galaxy sample used in weak cosmological lensing is 
selected by both magnitude and galax y size. Magnification due 
to len sing changes both quantities (e.g. Bart elmann & Schneider! 
1200 11) . and therefore introduces a correlation between number 
density and convergence. If the number density of galaxies with 
fluxes higher than some flux S , and sizes larger than some size 
R can be written as power laws, 



p(>S,>R)ozS'°R- 



(15) 



the observed number density differs from the unlensed one po to 
first order, becoming 



p = p (l + qic); q-Ka+p-l), 



(16) 



according to Schmi dFet alJ d2009allbh . This simple linear model 
assumes that the galaxy selection fu nction depend s on flu x and 
size in a separable way. We refer to IWvifhe et al.l (|2003) for a 
generalisation that takes into account correlations between the 
two quantities. 

Magnification and size bias induces a lowest-order correc- 
tion to the lensing power spectrum which has the same form as 
for reduced shear. Therefore, we can add the corresponding cor- 
rection term P^K to Eq. with 



P m ™(€) * AaYF mn (a k )Q M (€,a k ); 



k=\ 



Q fl (€,a) = qQ(£,a), 

where the function Q is given in Eq. (fT2l . 



(17) 



2.4. Fitting formulae 

For simplicity, we define Q (0) = Q and = dQ/dp a for 
a = 1 . . . M. These correspond to the ^-dependent terms in 
Eqs. ( fT2l and (1141 . which behave as near-power laws for both 
small and large (. With the abbreviation y = ln(, we perform 



two linear fits of In | Q^OOl for y < y\ ow = 3 and y > y up = 1 1 .5, 
respectively. We smoothly piece together these two asymptotic 
functions with a cubic spline s a (y) such that the composite func- 
tion In h^ a \y), a = . . . M, and its first derivative are continuous 



\nh (a) (y,a) = 



1 low 



(y, a) for y < y iow 



s^>(y, a) for y iov/ <y <y up 
8%\y, a) for y > y up 

for the linear functions 

8^(y,a) = b^(a)-y + bf\ay, 
g { ${y,a) = bf(a)-y + bf{a), 
and the cubic spline 



(18) 



(19) 



s (a \y,a) 




^-yiowj v . 


A 








= Si ( owCyiow,a) 




rf\a) 


= A-^(a); 




rfia) 


= 3gup ) ()'up, a) 


- A • b ( f(a) - 3r< a) (a) - 2rf ) (a); 


rf\a) 


= A • bf{d) - 


24W «) + ^\a) + r[ a) (a). (20) 



The ratios Q^'/h^ are then fitted by polynomials of order in 
y = ln£, at the fiducial cosmology p = po, and for each scale 
factor a = a k , k = 1 . . . N a , 

N 

Q (a \y, a) « h (a \y, a) ■ ^ c^ ) (a)y i ; a = . . . M. (21) 

;=o 

We note that we cannot fit the logarithm of the functions Q , 
since the derivatives of Q with respect to some cosmological 
parameters change sign. 

The fit coefficients b < f \cf f ^ are expected to smoothly vary 
with the scale factor a. We therefore perform another series of 
fits by polynomials of order Nb and A^ c , respectively 



7=0 j=0 



(22) 



The two matrices B (t?) and C ia) for a = . . . M defined in the 
last equation completely determine our approximation for the 
^-dependent terms of the reduced-shear power spectrum correc- 
tion. 

3. Results 

We use a fiducial flat ACDM cosmological model with 
WMAP7-like parameters, Q m = 0.2 7, fl b = 0.045, ft = 
0.71, o-g = 0.8 and n s = 0.96 dKomatsu et al.l l2010h . 
The dark-matter bispect r um i s calculated according to 
Scoccima rro & Couchmanl d2001l) . Their fitting formula is 
accurate to only 30% - 50% on small scales; this is however 
sufficient for our purpose. We note also that the bispectrum is 
not calibrated for any dark-energy model other than ACDM. 
For the m atter power spectr um, we use the 'halofit' fitting 
formula of Smith et al. (2003) and the transfer function 'shape 
fit' from Eisenstein & Hu (1998). Following the icosmo.org 
code (Ref regier et al.ll2008l) for models with w + — 1, we modify 
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'halofit' to interpolate between ACDM and w = -1/3, which 
beha ves in a similar way to an OCDM model (for more details, 
see lSchrabback et al]|2010l) . 

The function Q and the corresponding fit with the composite 
function given in Eq. OTT) are shown in Fig. Q] The polynomial 
that is part of the composite function has order N - 6. The fit- 
ting coefficients bf } and c' 0) are plotted in Fig. [2] Although the 
higher-order polynomial coefficients have relatively low ampli- 
tudes (right panel), we found that a polynomial of order 6 is nec- 
essary to provide a good fit to the reduced-shear power spectrum, 
as discussed below. For the polynomial fits of the coefficients as 
functions of a (l22l . we chose Nb = N c = 3. These cubic polyno- 
mials provide sufficient accuracy, in particular for a — > 1 , where 
the coefficients show the most variation. This is important be- 
cause the reduced-shear correction spectrum in Eq. fl% obtains 
a large contribution from large a. 

We perform the fits in the Grange between 0.1 and 2 x 10 5 . 
The functions Q (a ^ are not perfect power laws, therefore the fit 
for large t is not excellent. We found an improvement of our 
final results for the total power spectrum by adding 0.05 to 
after performing the fits. 

The fitting functions for the reduced-shear power spectrum 
corrections provide accurate results over a wide range in I. 
We illustrate the case of a single redshift bin with distribution 
p(z) = (z/zo) a ■exp[-(z/z ) /J ] and parameters a = 2,/3 = 1.5, and 
Zo = 0.5. The maximum considered redshift is z max = 2, which 
results in a mean redshift of 0.75. At I = 10 3 the reduced-shear 
correction to the convergence power spectrum starts to become 
important. On smaller scales, the latter underestimates the total 
power spectrum by more than 1%. For 10 3 < t < 10 5 , we fit 
Pg ' (Eq. to better than 20%. This is sufficient to provide an 
approximation of the total power spectrum at the percent-level. 

In Fig. [3] we the plot the ratio of the reduced-shear power 
spectrum derived using our fitting functions to that obtained by 
numerical integration of Eq. This is compared to the case of 
no correction for reduced-shear, corresponding to just the con- 
vergence power spectrum. In this latter power spectrum, a down- 
ward bias is evident, since the power is underestimated. This bias 
increases from 1% at I = 1000 to 6% at t = 10 s for the fiducial 
cosmology. In models with more structure, this bias is larger, 
e.g., 8% for erg = 0.93. In contrast, our fitted correction is accu- 
rate to better than 1% for £ — 2 x 10 5 at the fiducial point. 

We test different redshift distributions by changing to zo = 
0.3 and zo = 0.7, corresponding to mean redshifts of 0.45 and 1, 
respectively, and also a redshift bin of width 0.1 around z=0.75. 
In all cases, the fitting formula remains accurate to within 1 % at 
the fidicual model. 

In Table [1] we denote the parameter ranges for which our 
fitting formulae are accurate to better than 1% for I < 10 4 , and 
better than 2% for £ < 2 x 10 5 . 



4. Conclusions 

The lowest-order correction term for the reduced-shear power 
spectrum is the domina nt contribution from lensing-only effects 
dKrause & Hiratal [20091) . This term is proportional to the third 
power of the gravitational potential (f>, and includes an integral 
over the lensing bispectrum. In this paper, we have presented 
fitting formulae for this integral and its derivative with respect 
to cosmological parameters. This has permitted a more effi- 
cient calculation of the reduced-shear correction power spectrum 
within a parameter range comprising the 68% confidence region 
of WMAP7. The fitting scheme reduces the computational effort 
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Fig. 1. The function Q (0) (£) = Q{€) for different scale factors a as 
indicated in the panel (points), and the composite fitting function 
(Eq.EU lines). 

Table 1. Parameter limits for where the accuracy of the fitted 
reduced-shear power spectrum is better than 1% for £ < 10 4 , 
and better than 2% for £ < 2 x 10 5 , see also Fig. [3] Note that for 
each given range, the other parameters were kept fixed. 
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0.35 


2 


fide 


0.33 


1.03 
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W 


-1.6 


-0.6 
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Q b 


0.005 


0.085 


5 


h 


0.61 


1.11 


6 


OS 


0.65 


0.93 


7 




0.86 


1.16 



from tens of minutes for the full integration to a fraction of a 
second. 

For currently available weak-lensing data, the uncertainty in 
the non-linear power spectrum on small scales is still larger than 
the bias due to the reduc ed-shear. For exam ple, the widely-used 
fitting prescription by Smith et al. (2003) strongly underesti- 
mates the power on small scales. More recent numerical simula- 
tions however provide fit ting formulae that improve the accuracy 
by a f actor 5 to 10 (Heit mann et alj|201 Pi 120091 : lLawrence et al.l 
2010). Moreover, their emulation scheme can be applied to a 
wide range of cosmological parameters. This range is comprised 
in the hypercube for which our fits are valid. In combination, 
these results provide us with predictions for the weak lensing 
power spectrum that are accurate to the few percent level. 
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are polynomials in a or order N c = 3. Right: The power-law coefficients bf\ i = . . .3 of the composite function M 0) (Eq.[T8Tl. The 
lines are polynomials in a of order Nb = 3. 



& 10778725, 973 Program No. 2007CB 815402, Shanghai 
Science Foundations and Leading Academic Discipline Project 
of Shanghai Normal University (DZL805). 



B< 4 >: 



1.1983 -1.5973 3.3941 -2.4472 ^ 

21.942 28.174 -32.269 10.746 

-4.4314 9.5715 -14.092 6.4373 

U5.097 -124.48 200.67 -109.85 J 



Appendix A: Fitting matrices 

The matrices B {a) and C (Q ° d22|, a = 0...M contain the coef- 
ficients of the reduced-power spectrum fits given in Eq. (I2T1) . 
Here, we provide the numerical values from our fits. The index 
a = corresponds to the function Q fl% in the fiducial cosmol- 
ogy, and a - 1 ... 7 to its derivatives with respect to cosmolog- 
ical parameters (see Table Q]). The matrices are also available in 
electronic form with an example codtfl 



o(0) 



?(3) 



1.2157 -1.7061 3.613 -2.588 

22.054 14.525 -10.134 -1.2063 

-1.7221 -4.2894 6.9843 -3.7669 

Ul-719 34.891 -40.957 6.3019 

1.2963 -2.2435 4.7574 -3.4183 

24.802 15.438 -12.793 0.30538 
-0.42199 -14.769 28.968 -16.938 

62.035 -22.702 -10.113 15.12 

1.2314 -1.826 3.9322 -2.9291 

27.286 -8.0841 21.999 -18.426 

-3.2808 4.3912 -6.464 2.8446 

62.447 -73.471 122.62 -71.766 ) 

1.1772 -1.523 3.3478 -2.6269 ^ 

17.16 39.066 -56.363 24.476 

1.7985 -17.673 25.431 -12.429 

-1.2271 192.08 -259.8 113.95 



http://www2.iap.fr/users/kilbinge/reduced 



?(5) . 



q(6) 



?(7) _ 



C (0) = 



C 



■0) _ 



1.3643 -2.6937 5.6748 -4.0482 

22.681 17.078 -13.713 0.12838 

-1.6597 -3.5025 5.6229 -3.2007 

44.397 11.923 -3.5699 -10.617 



( 1.2149 -1.7012 

24.779 8.4417 

-2.5782 -0.084018 

53.073 -11.638 



3.6033 -2.5827 
-0.27702 -6.3675 
0.53439 -0.56922 
28.888 -27.294 



1.3469 -3.3556 6.6862 -4.5164 ^ 

25.294 6.7075 4.1411 -10.121 

-1.0794 -6.1604 9.2899 -4.7514 

37.736 49.648 -61.197 16.44 



0.56428 

0.12548 

0.1557 

-0.063989 

0.0087569 



2.3001 

-0.94677 

-0.65321 

0.43324 

-0.093102 



-3.9649 

2.59 

0.89795 

-1.0211 

0.24571 



2.427 

-2.04 

-0.58724 

0.7215 

-0.1693 



-0.00055709 0.008454 



-0.022712 0.015138 



1.4809 -10- 5 -0.00027429 0.00072501 -0.00046762 



-0.43501 -3.1459 5.6111 

-0.087793 0.80516 -2.761 

-0.10227 0.33267 -0.18646 

0.073105 -0.5126 1.2891 
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Fig. 3. The fitted reduced-shear power spectrum (thick lines) and the uncorrected convergence power spectrum (thin lines). Both 
quantities are plotted as ratios with respect to the reduced-shear power spectrum (|6), where the first-order correction Pi has been 
obtained by numerical integration of Eq. (TT~2b . In each panel, we show the power spectra at the fiducial model (dashed lines) and the 
lower (solid) and upper (dashed) boundaries within which a 2% accuracy is reached for € < 2 X 10 5 , see TableQ] 
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